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The Hamiltonian for a system of relativistic bodies interacting by their gravitational field is found in the post- 
Minkowskian approximation, including all terms linear in the gravitational constant. It is given in a surprisingly 
simple closed form as a function of canonical variables describing the bodies only. The field is eliminated by 
solving inhomogeneous wave equations, applying transverse-traceless projections, and using the Routh func- 
tional. By including all special relativistic effects our Hamiltonian extends the results described in classical 
textbooks of theoretical physics. As an application, the scattering of relativistic objects is considered. 
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Introduction. — The problem of motion of bodies under 
gravitational interaction has been the central issue in general 
relativity from its birth. As early as 1916, Einstein, Droste, 
Lorentz, and de Sitter started to develop post- Newtonian (PN) 
approximation methods. These are based on the weak-field 
limit in which the metric is close to the Minkowski metric 
and the assumption that the typical velocity v in a system di- 
vided by the speed of light c is very small, v/c ^ e. The 
deviation from the flat metric can be characterized by the 
Newtonian potential <I>; so for a binary system, for example, 
4>/c^ ^ v^/c^ ^ e^. In an appropriate limit (as e -^ 0), the PN 
approximation yields Newton's equations. In 1916 Einstein 
also worked out "the linearized approximation" to general rel- 
ativity in which the flat-space wave equations for the devia- 
tions hjjv from Minkowski metric with an energy-momentum 
tensor T^v of matter as the source are considered. This work 
marks the beginning of the post-Minkowskian (PM) approx- 
imation methods: the weakness of the gravitational field is 
assumed - the deviations of the metric from flat are small, but 
no assumption about slowness of motion is made. In an ap- 
propriate limit, when a suitable parameter, usually identified 
with Newton's gravitational constant G, is approaching zero, 
the PM approximation yields equations of special relativity. 
The "historical" period culminated in 1938 when Einstein, In- 
feld, and Hoffmann, by formal expansions in c^', investigated 
the n-th iterated field equations of the PN scheme. 

Since the 1950s numerous investigations of both PN and 
PM approximations appeared; see, for example, refined com- 
prehensive reviews |1, 2]. Because of the evidence that the 
orbit of the binary pulsar PSR 1913H-16 decays as a con- 
sequence of the emission of gravitational waves, the most 
promising candidates for the detectors such as LIGO, VIRGO, 
and GEO600 became binary neutron stars or black holes. This 
led to new studies of higher-order equations of motion. Exten- 
sive reviews summarizing these developments were prepared 
recently din. 

The Hamiltonian methods have been widely used in general 
relativity in such problems as perturbation of black holes, dy- 
namics of anisotropic cosmological models, etc. In the prob- 
lem of motion of gravitating systems the Hamiltonian yields 



equations of motion, expressions for the mechanical energy 
and angular momentum; in the case of binaries, the result 
for the last stable circular orbit follows from d^H/dr^ = 0. 
The Hamiltonian approach to the PN approximation was ini- 
tiated by Kimura and Toiya [5J and developed recently in the 
work of Schafer, Jaranowski and Damour (see j^ for refer- 
ences). The canonical formalism of Arnowitt, Deser, and Mis- 
ner (ADM) |6] is commonly used. 

Before we enter into details we wish to indicate why our 
final result-the closed form Hamiltonian (fTTI ) including all 
terms linear in G — is of importance in the theory of gravity 
and can be relevant in other branches of theoretical physics 
and astrophysics: (i) Since PM approximation can be re- 
stricted to slow motions, the Hamiltonian describes PN ap- 
proximations to any order in 1/c when terms linear in G are 
considered. When considering only two particles and terms 
up to the second order (< v^/c^) in center-of-mass system we 
precisely recover the Hamiltonian given by Landau and Lif- 
shitz [7], except for the last term which is ^ G^. Moreover, 
we checked that, after a suitable canonical transformation, it 
yields precisely all the terms linear in G in 3PN approxima- 
tion of [8]; and we calculated corresponding terms in 4PN 
order which have not been given thus far. These will be pub- 
lished elsewhere, (ii) An electromagnetic counterpart of our 
Hamiltonian was derived by completely different methods by 
Kennedy |9]. The first terms in the expansion in c^^ of the 
Kennedy Hamiltonian correspond to the well-known Darwin 
Lagrangian (see e.g. [7, 10]) but Kennedy's treatment goes 
beyond Darwin's approximation. Our procedure shows how 
the Darwin Hamiltonian or Lagrangian can be generalized in 
both gravitational and electromagnetic cases. The Darwin 
Lagrangian or Hamiltonian have many applications in both 
classical and quantum domains (see, e.g., ||10l]). The extension 
for two point charges considered up to order c^^ (included) 
has been given by the Golubenkov-Smorodinsky Lagrangian, 
which is described in |7] and discussed in jllll . (iii) As our 
Hamiltonian can describe particles with ultrarelativistic veloc- 
ities or with zero rest mass, it is well suited for treating grav- 
itational scattering of such objects. These are closely related 
to gravitational scattering of shock waves. At high energies 



in quantum scattering processes, the interaction is well ap- 
proximated by the classical collision at the speed of light, as 
first shown by the influential work by 't Hooft (see 1 12] for de- 
tailed description and methods), (iv) Because of its simplicity, 
our Hamiltonian yields a convenient starting point in investi- 
gations of gravitational lensing by extended gravitationally in- 
teracting objects since it provides a unified description of both 
gravitating bodies and deflected photons or neutrinos, (v) Our 
method of deriving the Hamiltonian by using the Routhian 
and the induced canonical transformation, etc. should be of 
interest in other branches of theoretical physics. 

The Hamiltonian approach to the PM approximation was 
first undertaken in 1986 [13]. Following the ADM canoni- 
cal formalism, the gravitational field is described by hj-^ , the 
transverse-traceless part of /i,y — gjj — 5,y (hj^^ = 0, hfj'^j — 0, 
i,j = 1,2,3), and by conjugate momenta Jt'-'^^ . The system 
of bodies located at x„, a = l,...N, with rest masses nia and 
momenta pa, has the energy and linear momentum densities 

73 T^'^n^ny = I„ {g'JpaiPaj + ml) ^ 5(x - x„), ~7^- T, ^n^ = 
'LaPai5{x - Xa), where 7 = det{gij), g'J is inverse to gij, n^ is 
a unit timelike normal to hypersurface x** = const, and T^^ is 
the energy-momentum tensor of the matter system. Hereafter, 
we call its constituents the "particles", but they may well rep- 
resent neutron stars or black holes. This is substantiated by 
"general relativity's adherence to the strong equivalence prin- 
ciple": black holes and other bodies obey the same laws of 
motion as test bodies; see, e.g., yj 



It is convenient to choose four coordinate conditions n" = 
and gij = (1 + ^0) 5ij + hff. The standard ADM Hamil- 
tonian (cf. [[^), H = {l/l67zG)§dSi{gij,i~gjj,i), then be- 
comes, using the Gauss theorem, H — -~{l/167zG) fd^xA<j). 
The integrand A(j) = djdiCJ) can be expressed in terms of x^, p„, 
hjj^ and n'^^^ from the constraint equations. By expansions 
in powers of G and after adopting suitable regularization pro- 
cedures of integrals involved (see the Appendix in fO*]), one 
can determine the Hamiltonian. The Hamiltonian worked out 
in 111 311 includes terms ^ G^. 



Here we start from the same Hamiltonian neglecting, how- 
ever, terms ^ G^. Still, the form of the Hamiltonian is quite 
complicated [see Eq. ([U]: rather than just quantities asso- 
ciated with particles it involves field variables, non-local TT 
projections, and integrals. The main purpose of this Letter 
is to show that due to somewhat magical simplifications the 
fields hjj^ , assuming they are generated just by particles, can 
be expressed entirely in terms of particles' variables, and after 
two canonical transformations, the Hamiltonian can be cast 
into the closed form (fTTT l. It is local in particles' canonical 
variables, involving no field variables. It is exact up to terms 
linear in G, and within the same accuracy it yields equations 
of motion for particles moving with possibly ultrarelativistic 
speeds, including those with zero rest mass. 

1PM Hamiltonian. — The Hamiltonian describing A^ parti- 
cles (Xa, Pa) and their gravitation field (/i/J-^, n'-'^^) accurate 



up to the terms linear in G reads 

Min=VOT„--G> l + =2+=2 (1) 

+ tG 2^ — (7p„.pi + (p„.n„/,)(pi.n„i)) 
-\V-^hJl(x = x..) 



where ot„ = (m^ + p^)', na^r,,/, = x„ -Xf,, |n„/,| = 1. [This 
is Eq. (18) in 111311 . However, although we also put c = 1, 
we keep explicitly G and do not put 167rG = 1; the "geomet- 
ric" momenta tt'-'^^ are thus different from the true Tt'-i^^ ~ 
{l6nG)^^K'^^^ . Notice that the regularization is also needed 
in the terms containing hj-'^ , n'^^^ , and hj-^ {x ~ x„).] The 
equations of motion for particles are standard Hamilton equa- 
tions. The Hamilton equations for the field are 
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(2) 



here the variational derivatives and the TT-projection operator 

^ki^''' = 2 {AikAji + AiiAjk~AijAki)A^^, Atj = 5ijA~didj, 
appear. These equations imply the equations for gravitational 
field in the first PM approximation to be the following wave 
equations: 
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16;rG5/^™y:^^5(3)(x-x„) 



(3) 



As is usual in the system of particles and field, the field is 
source of particles' accelerations and particles are sources of 
the field. Once O(G^) terms can be neglected, however, this 
coupling simplifies. Since both the field and the accelerations 
Pa are proportional to G, the changes of the field due to the 
accelerations of particles are of the order 0{G^). Therefore, 
we can assume the field to be generated by unaccelerated par- 
ticles; it is thus given only by their instantaneous positions 
and velocities. If such a field is used in equations of motion of 
particles, the gravitational interaction becomes an "action at a 
distance". 

Given the linearity of the field equations (O, we can write 
hf/{'>^) == Lfo/^^-^(x;Xfo,pfo,Xfo), where hjf {x;xt,pt,'>ib) is the 
contribution to the field at point x generated by the particle 
moving at x^ with velocity x;,. Because the d'Alembertian 
operator commutes with the TT-projection operator, the field 
from each particle involves the solution of the scalar wave 
equation with point source - the boosted static spherical field 
^ 1/r. Denoting x — x„ = n„|x — x„| and cos0„ = n„.x„/|x„|, 
the solution of (O can thus be written as 



/.^^(x)=^"'i: 



AG 1 



PhkPbl 



nib |x-Xf,| 



1 — x^ sin^ db 



(4) 



The action of 5jJi in equation (|4|l consists of two steps: 
first, one has to solve the Poisson equation twice and then 
evaluate a number of partial derivatives. The first step is fea- 
sible due to the form of the boosted spherical potential. It can 
be shown to satisfy the relation (here v ~ |xa|) 



2\2 



3(1 -v/) 



l-v2sin^a, 



7 ,L .2 ^^ 
dv dv^ 



(5) 



|x-x«| v'l-v2sin2 0„ 



Hence, instead of working out A '^, we can evaluate 
[8 + 7v d/dv + v^ c/2/(iv2j ^ j g., the elliptic partial differen- 
tial equation of the fourth order can be simplified into an in- 
homogeneous linear second order ordinary differential equa- 
tion. It has a unique solution which is regular at v = 0. Even 
though this solution is quite complicated and contains loga- 
rithmic terms such as In cos da, they cancel out when partial 
derivatives are combined. After somewhat lengthy calcula- 
tions (the details of which will be given elsewhere), we find 
the field of a moving source 



^7^(x;Xi,Pi,Xi) 



G 

|x-x;,| nib 



^yOT^i bPfo- (n/'-Pfc)'(33' + 2)] 5,7 + 2 [l -x2(l -2cos2 0,)] p,ip,j 

[(2 +y) (ni,.pfo)2- (2 + 3^ - 2x^) pf ] nbitihj + 2(nfo.pfo) (l - x^ + 2y) {tibiPbj + Pbifibj) } 
0{mbiib~9b)G+0{G^)- 

I 



(6) 



here y = jb 



1 - x2 sin2 db. As indicated by the symbol tained. The field part of the Routhian 



Oirfibiib ~ Pb), the last expression gets simplified by using 
Xft = Pb/^b and we anticipate that later Oimbii-h ^ 9b)G will 
turn into terms ^ 0[G^). 

In order to later suppress field degrees of freedom, we shall 
turn to the Routh functional (see, e.g., ISD) 



R{lLa,VaMjMj)=H 



1 



\6nG 



d^^TT.Ji^TT 



(7) 



which is "the Hamiltonian for the particles but the Lagrangian 
for the field." While the functional derivatives of the Hamilto- 
nian yield the time derivatives of the canonically conjugated 
field, the functional derivatives of Routhian vanish if the field 
equations ^ hold. Their solution (|4]i is non-radiative and can 
thus be substituted into the Routh functional without changing 
the Hamilton equations for the particles. Hereafter the symbol 
hjj'^ is a shortcut for the solution of (O depending on coordi- 
nates, momenta and velocities of the particles. So the reduced 
Routhian (7) which becomes Hamiltonian i/(xa,p„,x„) is ob- 



Rf 



1 



leKG 



d'-7ih[^h[l 



,TT 



-h^firfl) 



(8) 



still needs to be transformed into an explicit function of the 
particles' variables. Using Gauss's law in the first term and 
integrating by parts the second term, we arrive at 



Rf = -■ 
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(9) 



64nG dt 



d'xih'JhU) 



The field equations ^ imply that the first integral (in which 
the self-interaction term is discarded) directly combines with 
the "interaction" term containing Y-^na^ Pai Paj hJF i'^a), so 
only its coefficient is changed. Remaining terms do not mod- 
ify the dynamics of the system in our approximation. The 
Hamiltonian thus takes the form 
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-^ — hij (x = Xfl;Xft,pi„Xi, . 



jG Y, — (7pa.p/, + (p„.n„i)(pfo.n„i)) 



(10) 



Dropping out the total time derivatives in (|9]l means a canonical transformation, but the new canonical coordinates will not 
be denoted by primes. In fact, another change of coordinates has to follow since the Hamiltonian ( fTOl i is a function of x,,. 
We define new momenta by putting p',; = pai — dH/dxai, and then eliminate Xa by introducing new Hamiltonian //'(xa,pj,) = 
//(xa,pa(p',),x„(p^)) — LiJCi)i(Pa) dH/dxbi- Since dH/dxa ^ G, the only change in the Hamiltonian which is linear in G comes 
into the kinetic term ^m^ from the last change of momenta. This change is exactly cancelled by the sum Y.h^bi{l?a) dH/dxbt in 



i/'(xa,p^). We now make simple substitutions i^ 



dm„ 



f , COS 0„ ^r^,yl = \ -x2sin2 0, ^ 1 +m„-2 



(na.Pa 



and, again, omit primes. At this moment we can substitute for hj^ in ( fTOl i the solution (|6]l, in which the above substitutions turn 



^2ll 



the term 0(m/,Xf, — pi)Ginto term ~ (9(G). In this way, using the shortcut y/^a = niy [m^ + (nt„.pi) ]2, we finally arrive at the 
Hamiltonian in the first post-Minkowskian approximation: 
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(11) 



2(2(p„.pi,)2(p^.nfc„)2-2(p„.n/,„)(pfe.nfea)(p«.Pfo)pi + (p«.ni„)^p^-(p„.p,,)2p^)=2 
2 [-p^(Pi.n;,„)2 + (p„.nfoa)2(pfo.ni„)2 +2(p„.nfoa)(pfe.nfoa)(p„.pfo) + (p„.pi)^ - {9a-nbafVb\ 
[-3p,^(Pfo.ni,fl)^ + (p„.nfo„)^(Pfo-nfe,)^ + 8(p„.nfca)(p/,.ni„)(p„.pfo) + p2p^-3(p„.nfe„)2p2]3;^„ 



This is our main result. The Hamiltonian for a many-particle 
gravitating system in post-Minkowskian approximation, i.e., 
including all terms linear in G, was derived in the closed form 
entirely in terms of the variables of the particles. Putting 
G = it becomes standard Hamiltonian of A^ noninteracting 
particles in special relativity. 

Let us yet note that in [ 14] the post-Minkowskian action for 
a helically symmetric binary solution was constructed; how- 
ever, it turns out that due to the restriction to helical worldlines 
ambiguities may arise. We will return to this issue elsewhere. 
It should be also useful to derive the Hamiltonian (fTTT l within 
the effective field theory approach to gravity 111511 . 

Scattering. — As an application of the Hamiltonian obtained 
we calculate gravitational scattering of two possibly ultrarela- 
tivistic or zero-rest-mass particles. When only terms linear in 
G are considered, the transferred momentum can be computed 
as Api — X^^Tpi dt , integrating along the straight line trajec- 
tories of noninteracting particles [pi is determined from the 
Hamilton equations using (fTTT il. If perpendicular separation b 
of trajectories (|b| is the impact factor) in center-of-mass sys- 
tem (pi = — p2 = p) is used, p.b = 0, we find, after evaluating 
a few simple integrals, that the exchanged momentum in the 
system is given by 



Ap 



b G 



— 2 — 2 
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The quartic term is all that remains from the field part hj^ . 
It is not difficult to show that ( fT2] i agrees with the result 11 1611 
obtained by a very different method. 

The authors benefitted from the exchange program between 
Charles University, Prague, and Friedrich Schiller University, 
Jena. TL. and J.B. acknowledge the partial support from 



SFB/TR7 in Jena, from the Grant GACR 202/06/0041 of the 
Czech Republic, and of Grants No LC 06014 and the MSM 
0021620860 of Ministry of Education. J.B. is also grateful for 
the support of the Alexander von Humboldt Foundation. 



[1] T. Damour, in Gravitational Radiation, edited by N. Deruelle 
and T. Piran (North-Holland, Amsterdam, 1983), pp. 59-144; in 
Three Hundred Years of Gravitation, edited by S. Hawking and 
W. Israel, (Cambridge University Press, Cambridge, England, 
1987), pp. 128-198. 
[2] K. Thome, in Three Hundred Years of Gravitation (Ref. [1]), 

pp. 330-458. 
[3] L. Blanchet, Living Reviews in Relativity 9, 4 (2006) 

http : // www.livingreviews ■org/lrr-2006-4, 
[4] T. Futamase and Y. Itoh, Living Reviews in Relativity 10, 2 

(2007), http : /Twww . livingreviews . org7Trr-2007-2j 
[5] T. Kimura and T. Toiya, Prog. Theor. Phys. 48, 316 (1972). 
[6] R. Amowitt, S. Deser, and C. W. Misner, in Gravitation: An 
Introduction to Current Research, edited by L. Witten (Wiley, 
New York, 1962), pp. 227-265. 
[7] L. D. Landau and E. M. Lifshitz, The Classical Theory of 

Fields, 4th ed. (Pergamon, New York, 1975). 
[8] T. Damour, P. Jaranowski, and G. Schafer, Phys. Rev. D 62, 
021501(R) (2000); Phys. Rev. D 63 029903(E) (2000); Phys. 
Lett. B 513, 147 (2001). 
[9] F. J. Kennedy, J. Math. Phys. (N.Y) 16, 1844 (1975). 
[10] J. D. Jackson, Classical Electrodynamics, (John Wiley, New 

York, 1975), 2nd ed. 
[11] T. Damour and G. Schafer, J. Math. Phys. (N.Y.) 32, 127 (1991). 
[12] P. D. D'Eath, Black Holes: Gravitational Interaction, (Claren- 
don Press, Oxford, England, 1996). 
[13] G. Schafer, Gen. Relativ. Gravit. 18, 255 (1986). 
[14] J. L. Friedman and K. Uryu, Phys. Rev. D 73, 104039 (2006). 
[15] W. D. Goldberger and I. Z. Rothstein, Phys. Rev. D 73, 104029 

(2006); Gen. Relativ. Gravit. 38, 1537 (2006). 
[16] K. Westpfahl, Fortschr. Physik33, 417 (1985). 



